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1 Introduction 

It has been proved that on 2- dimensional orientable compact manifolds of genus 
g > 1 there is no integrable geodesic flow with an integral polynomial in momenta, 
see m |[. There is a conjecture that all integrable geodesic flows on T 2 possess an 
integral quadratic in momenta (for more details see 0). All geodesic flows on S 2 
and T 2 possessing integrals linear and quadratic in momenta have been described in 
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]5|, |1J, using some ideas going back to Darboux ]3[]. So far there has been known only 
one example of conservative system on S 2 possessing an integral cubic in momenta: 
the case of Goryachev-Chaplygin in the dynamics of a rigid body, see ||. 
The aim of this paper is to propose a new one-parameter family of examples of 
complete integrable conservative systems on S 2 possessing an integral cubic in mo- 
menta. We show that our family does not include the case of Goryachev-Chaplygin. 
The paper is organized as follows: 

We start with the investigation of the following initial value problem 
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-7u"u'r 3 - 2u" 2 r 4 + u"ur 2 - 2u' 2 r 2 + u'ur + u 2 , 



tt(l) = 0,« , (l)=l,M /, (l)=T-l, (1) 

where u : R + h R, r h u(r), and prove that there exists a positive constant T 
such that the one-parameter family {\l> T } re ( T ) of solutions of this problem defines 
a one-parameter family of smooth conservative systems on S 2 with energy 

H T (r, <p t dr, dip) = lA %^) ~ ^ + Kr ~ *') ^ cos * (2) 

in polar coordinates r and ip, possessing nontrivial additional integrals cubic in 
momenta. 

We note that the limit point of our family for r = corresponds to the function 
\?o = \ (r — -) which defines by (0) the standard metric of constant curvature 
K = 1 on S 2 . " 
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2 Existence of C°° solutions 

In this chapter the problem of the existence of the smooth solutions of ([!]) on (0, +00) 
will be considered. 

Initial value problem ([!]) can be replaced by the initial value problem 

x'x'" = xx" - 2x" 2 + x' 2 + x 2 , x(0) = 0, x'(0) = 1, x"(0) = r. (3) 

This reduction is accomplished by the change of variable r — > t — logr which 
changes ([J) to (^). Initial value problem (||) is autonomous and more convenient for 
further investigations. For this reason, the theorem which follows will be stated for 
(H), but proved for (|3]). 

Theorem 2.1 There is a positive constant T such that solutions of (|7|) exist on 
(0, +00) for r E (— T, T) . 

Proof. Initial value problem (f|) has a unique solution x(t) = O r (t) which is positive 
on (0,e) and negative on (—£,0) for a sufficiently small e. 

Let us consider the case t > 0. 

The differential equation from ([|) can be replaced by the following differential equa- 
tion of the second order 

q = - (1 + 2q 2 - 3q 4 + q - 7q 2 q - 2q 2 ) 
q ' 

with function q(t) = R'(t) where R(t) = \ogx(t). We may rewrite this equation as 
a system of differential equations of the first order: 

q = p, p=-(l + 2q 2 -3q 4 +p- 7q 2 p - 2p 2 ) . (4) 
q v 

Since system (f|) is symmetric with respect to q 1— > — q, 1 1— > —t, it suffices to consider 
the case q > 0. 

In order to obtain the phase portrait of (£|) we may consider the following smooth 
system 

q = qp, p = 1 + 2q — 3q + p — 7q p — 2p . (5) 

The solutions of (|4]) are obtained from the solutions of (|5|) by a reparametrization. 
System (||) has four singular points: two saddle points p = l,q = Q, p = —~,q = 
and two knots p = 0, q = ±1. 

The aim of our further investigations is to show that the orbits of @, corre- 
sponding to the solutions of (J3|), when r belongs to a certain interval, converge to 
the singular point q = l,p = 0. 

For any solution x{t) = T (t), t G (0, e) of (||) there is an orbit {T T : p = p T (q)} 
in the phase spase of (f|). 

For any solution of (|3]) it holds that R(t) — > — 00 as t —>■ 0+ and, therefore, q(t) = 
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R'(t) — > +00, p(t) = q'(t) — > —00 as t — > 0+. Thus p T (?) - ► _ 00 as g — > +00. 
We may now consider only the orbits of (|j) where p — > —00 as g — > +00. 
Show that if T\ > r 2 then p T1 {q) > p T2 ( < l)- Indeed, for a solution x(t) = 6 T1 (t) of (|3]) 
it holds 

Note that the function x(t) = sinht satisfies (|j) when r = 0. The related orbit of 
(|) has then the following form {T : p = po(q) = — g 2 + 1}. 

So, the orbits of (f|), corresponding to the solutions of (|]), for r > converge to the 
singular point q = l,p = 0. 

We prove below that orbit (*), where p — > — | as g — ► 0+, corresponds to the 
solution of (Q) when r is equal to a negative constant — T and all orbits of (f|), lying 
between (*) and r correspond to the solutions of (§) for r 6 (— T, 0). 

Assume first that there exists a constant to < such that orbit r ro does not 
converge to the point q — l,p — 0. 

Consider the set W of the orbits of ([|), lying between r ro and T . Show that for 
any orbit in W the value q becomes infinite in a finite time interval. 
In fact, for any solution of (W) holds 

q(t) dq 

— = t + const. (7) 
p 

For orbits in W we have p < —q 2 + 1. Hence, the left hand side of (0) is bounded 
for any orbit in W as q —>■ +00. Without loss of generality we assume that for these 
solutions of (f|) t vanishes as q becomes infinite. 
We conclude that for any orbit from W it holds by @ 

{r + o{l))q < p + q 2 < 1 

and for a corresponding solution of (||) we get for t —> 0+ 

/ x X "(t) 

Thus, for a corresponding solution of (|3|) the function (logx'(t))' is bounded in an 
interval, containing 0+. Hence, lirn t _ +0 + x'(t) is bounded for a solution x(t) of (||), 
corresponding to an orbit in W and, therefore, x(t) = — > as £ — » 0+. 
Let us consider orbit (*) where p = p*{q). If (*) is the same as r T0 then T = — r . 
If (*) is not the same as r T0 then it belongs to W and, hence, corresponds to the 
solution of @ for r = —T, where T is equal to a positive constant and, moreover, 

T < -T . 

Thus we have shown that for any r G (— T, T) where T is a positive constant sytem 
(0) has a solution for all t > 0. 
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Let us assume now that there is no such orbit r ro and, so, all orbits of (J|), 
corresponding to the solutions of (j|), converge to the singular point q = l,p = 0. 
Clearly, orbit (*) corresponds to a solution x(t) = r)(t) of the differential equation in 
(§). As mentioned above, the function q(t) = ^M- becomes infinite in a finite time 
interval of t, say, as t — > 0+. Since p*(q) < —q 2 — kq for any k and for sufficiently 
large q, for r](t) it holds 

li m = _oo. (8) 

t^0+ 7]'(t) V ' 

Note that 77(0) is finite because R(0) - R(t) = J t °q(t)dt < for t > and 77(f) = 
expi?(t). So, there are two cases: 77(0) = or 77(0) 7^ 0. 

Assume that 7/(0) 7^ 0. It follows that t/(0+) = +00 and from (|) we get ?7"(0+) = 
— 00. 

Rewrite now the differential equation in (|3|) in the following form 

— = - 3 — - 2 ~ i~ + -• 

nn ryi rprpl rp 

Therefore, it holds 

rf = 3 P*(g) + q 2 - 1 _ 2 (p*(g) + g 2 -l) 2 + 
77 (7 q 

We obtain ?7'"(0+) = —00. This is a contradiction 7/"(0+) = —00. 

Assume that 77(0) = and 77^0) 7^ 0. From (|j) we get ?7"(0+) = —00. Rewriting 

the differential equation in (^) in the following form 

(rpM rp \ rp ( rfJ^ rp \ 2 

x' x' 

we obtain again ?7'"(0+) = — oo. 

Therefore, we only have to consider one case 77(0) = and r]'(0) = 0. Taking into 
account that rf{t) = rj(t)q(t) and q(t) > as t — > 0+, we conclude that rf(t) — > 0+ 
as t — > 0+ but on the other hand from (|8[) it follows that r]"(t) < 0. 
These contradictions finally show that there is an orbit r ro which does not converge 
to the point q = l,p = and, therefore, for any r G (— T, T), where 

p*(q) + g 2 
T = — hm < 00, 

solutions of (^) exist on [0, +00). 

Since a solution x(t) = Q T (t) of (§) equals — 0_ T (— £) if t < 0, solutions of (§) exist 
on (—00, +00). □ 



3 Asymptotic behaviour at infinity 

In this chapter we deal with the properties of the solutions of ([D for r G (— T, T) 
at r = and r = +00. 
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In this chapter initial value problem (|1]) is oft replaced by the autonomous initial 
value problem (^) and another initial value problem which is more convenient when 
r — > or r — ► +00. 

Theorem 3.1 There is a constant T such that for any solution \1/ T (r) of (|I[) with 
r G (— T, T) it holds 

+ %r - * T ) = /i r (r 2 )r = ^(^)-, 

where the functions £ T , ( T , v T , fi T are of class C°° and £ T (0) 7^ 0, Cr(0) 7^ 0. 

Proof. Initial value problem (|l|) can be replaced by another initial value problem. 
The change of variables r — > s = r -2 , it — > (7 = r~ x w changes (0) to the following 
initial value problem 

</" = 3^|V + S ^V> fC 1 ) = >4'(1) = -^"(1) = r (9) 
9 - 5 - 2g's 2 

We will prove that solutions of @ for parameter r G (—T,T) satisfy the regular 
initial conditions at s = 0: g(0) 7^ 0,</(0) < 00, < g"(0) < 00. 

Compute now for a solution of (|3|): 

x'{t) = {e W t){g-2g's), (10) 

x"(t) - ar(t) = 4exp(-3t)g"(s), (11) 
remember that s = exp (— 2t). 

Let us consider system (f|). Since the eigenvalues of the corresponding linear system 
are equal to —2 and —4, there exist functions P and Q of class C l such that 

q = Q(s) = 1 + Cs + o^s), 

p = P{s) = -2Cs + o 2 {s), 

where C is a constant, see ||. 
We may write: 

R{t) = J q = t + i; 1 {s). 

Let us show that ipi is of class C . By differentiation we obtain 

djhjs) = q - 1 = Cs + o^s) c0 
ds 2s 2s 

Write now for a solution of 

x(t) = exp R(t) = (exp t) exp ipi(s) = (expt)g(s), 

where we have used that g = r~ x u = exp(—t)u = exp(—t)x(t), and, therefore, 
<7(0) = exp Vi(0)^0. 
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For a solution of (|3|) which can be extended into infinity it holds 

x'(t) = x(t)q(t) = (expt)(expV>i(s))(l + Cs + 01(a)) = 

= (expt)g(s)(l + Cs + o 1 (s)), s = exp(-2t). 
On the other hand, from ( |IUD we obtain 

g(s)(l + Cs + o 1 (s)) = g(s)-2sg'(s). 

Thus, for any solution of ([|) where the corresponding orbit of converges to the 
singular point q — l,p — it holds g'(0) < oo. 

Let us show that for these solutions it holds < g"{0) < oo. 
We compute 

x"(t)-x(t) = x(t)(p + q 2 - 1) = (expt)(exp^i(s))(P(s) + Q 2 (s) - 1) = 
= (expt)(exp^(5))(-2Cs + o 2 (s) + (1 + Cs + 0l (s)) 2 - 1). 

Thus, 

= (expt)^(s), 

where ^2 G C 1 , ^2( s ) = °3( s )- Rewrite the differential equation in ([3]) in the 
following form 

I ^ 3/ — O ( X j X Zi\X X ) 

We obtain either x"(t) — x(t) = or 
(log(ar"(t) - x(t)))' = -3- - 2(x"(t) - x(t))(exptexpV»i(s)(l + Cs + oi(s)))" 1 = 

= -3- - 2^ 2 (s)(exp Vi(s)(l + Cs + oi(s))) -1 . 
Q 

So, 

(log(x"(t) - x(t)))' = -3(1 - Cs) + o 4 (s). 

Then it follows 

(log ((expt)^(s)))' = -3(1 - Cs) + o 4 (a). 

Thus, 

(< + log^(s))' = -3(1 - Cs) + o 4 (s). 

By integrating we get 

3C 

(logV^(s)) = -4t - — s + ^s(s), 

where ^3 G C 1 . 
Thus, 

3C 

ip 2 {s) = (exp(-4t))exp( — — s + tp 3 (s)) = s 2 tp A (s), 

since s = exp(— 2t). So, -0 4 G C 1 and ^4(0) = const 7^ 0. Taking into account (pT|) , 
we obtain g"(0) = ^4(0) and therefore < p"(0) < 00. 

So, we conclude that the solutions of (0) for r G (— T, T) are of class C°° in zero. 
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Let u{r) = \I/ r (r) and x(t) = Q T (t) be solutions of ([I]) and (^) respectively. Then 
we get that ^ T (r) = r (logr) and 

^(r) = e'(logr)i = - 2^(1)1 = C.(^), 

where ( T is of class C°°. From (|T0|) and the condition <?(0) 7^ it follows that 

Cr(0) ^ 0. 

Using ([!(]) and (pTTJ ) we may compute 

(^y + ^> - * T ) *; 2 r 2 = (e^(t) - e T (t)) e; 2 (t) = 

4exp(-t)/( S )(^ - 2^'( S ) S ) 2 = 4- ^(1) - 2^(1)1) g"{-) = u T (\)-, 
where v T is of class C°°. 

Since a solution x{t) = Q T (t) of (||) equals — 6_ T (— t) if t < 0, in the same way, we 
obtain 

= ^.(r 2 ), 

where £ T G C°°, f T (0) 7^ and 

(^> 2 + *J.r - tf r ) ^; 2 r 2 = £t T (r 2 )r, 
where /i T is of class C°°. □ 

4 Smooth conservative systems on S 2 

In this chapter we propose new examples of complete integrable conservative systems 
on S 2 in terms of solutions of (JlJ) . 

Theorem 4.1 The one-parameter family of conservative systems with energy (fjj, 
where \1/ G {\I/ r } Te (o,T) <™<^ a solution of (^), is a one-parameter family of 

smooth systems on S 2 possessing an integral cubic in momenta which is nontrivial, 
i. e. there is no quadratic or linear integral. 

These systems cannot be obtained one from another by a change of variables. 

This family of integrable conservative systems on S 2 does not include the case of 

Goryachev-Chaplygin. 

Proof. Prove now that the kinetic energy of (|2|) for any r G (0, T) is a C°° metric 
on S 2 . Using theorem |3.1| , we may write the kinetic energy as 

r 2 d<p 2 + dr 2 f 2 d(p 2 + df 2 
Ur 2 ) = CAr 2 ) ' 

where f = -, (p — —(p. Since functions £ T and ( T are of class C°° and, moreover, 
£ T (0) 7^ 0, Cr(0) 7^ 0, it follows that the kinetic energy is a smooth metric on S 2 . 
The potential of (§) also can be written in polar coordinates: 

V T = fi T (r 2 )r cos (p = u T (f 2 )r cos (p. 
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Thus, the potential is also a smooth function in r cos tp, r sin (p, < r < oo and in 
fcos<p, f sin<£>, < f < oo and therefore, a smooth function on S 2 . So, any system 
with energy (0) is a C°° conservative system on S 2 . 

We prove now the integrability of these systems. 

As mentioned above, if \I/ T satisfies ([I]) then the function 6 T = \I/ T oexp satisfies (H). 
In the coordinate system y = logr, (p = tp energy has the following form: 

h t = ®r ,2 {y){pl + P 2 y ) - e T ' 2 (y)(e T "(y) - e T (y)) cos ^, (12) 

where 6 r G {e T } re ( ,T)- 

We will prove that the following polynomial cubic in momenta 
3 

F r = P% + g (®Ay) cospp v - Q' T (y) sinppy) = p^ + E T (p,y,p v ,p y ) 

is an integral of the system with the Hamiltonian (0), i.e. {F T , H T } = 0. 
Denote by H T = Q T ' 2 (y)(p 2 + P 2 ) the Hamiltonian of the geodesic flow of metric 

- 2 1 

dSr = e /2 (y) ^ + dy ^ ' 

Write now energy as H T = H T + V T where the potential V T has the following form 

V T = -Q T ' 2 (y)(Q T "(y) - Q T (y)) cos^. (13) 

Note that {H,p^} = 0. In order to prove that {F T ,H T } = we must prove that 

{K,Py + {H T ,E T ((p,y,p v ,p v )} = and {V T , E T ((p, y,P(p,p y )} = 0. 
Compute 

{Vt,pI} + {Hr^E^ip.y.p^py)} = -36 T /2 (6 T " - 6 r ) siny?p 2 - 36 r ' 2 (6 r sin (pp 9 + 
6/ cos 'ppy)p ip + 36 T /2 (0 r ' cos (pp 9 - <d T " sin <pp y )p y + 30 T /2 T " sin y?(p 2 + p 2 ) = 0. 

By a computation we obtain {V T , E T (ip, y,Pip,p y )} = 




<p sin pQ' 2 (q' t q™ - e 2 - e^e r - e; 2 + 2ef ) = o. 



Let us assume that a system of this family has an integral which is independent 
of the energy and which is a polynomial of second degree with respect to momenta 
(clearly, this assumption includes the case of linear integrals). So, there is an integral 
F T of ( |12D which is quadratic in momenta. Thus, F T = A T {j> v ,p y ,tp,y) + B T (p,y) 
where A T {p v ,p y ,p,y) is a polynomial of second degree with respect to momenta. 
We may write {F T , H T } = {Ar^p^jPy, tp, y) + B T (<p, y), H T + V T } = and, therefore, 

{Ar(P<p,Py,(p,y),H T } = 0. 

- 2 

Thus, the geodesic flow of ds T has an integral which is a polynomial of second 
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degree with respect to momenta and, therefore, in some Liouville coordinate system 
u, v it holds 

A, = * ,M + vll v, = + (14) 

oti(u) + a 2 {v) a^u) + a 2 {v) 

where «i, a>2, 03, «4 are some functions. 

Geodesic flows on S 2 possesing an integral quadratic in momenta have been de- 
scribed in ||. In particular, it has been shown that Liouville coordinate systems 
for all metrics on S 2 of nonconstant curvature are unique up to linear combinations. 
Noting that (|T3| ) for any r G (0, T) has not form ( |HD in any coordinate systems 
ay + b, cup + d, where a, b, d are some constants, we prove that no system from our 
series possesses an integral quadratic in momenta. Moreover, since the solutions 
of (|3]) for r G (0, T) cannot be obtained one from another by a linear change of 
variables, systems with the energy (0) also cannot be obtained one from another by 
a change of variables. 

Prove that our family of smooth conservative systems on S 2 does not include the 
well-known case of Goryachev-Chaplygin from the dynamics of a rigid body. 
In the case of Goryachev-Chaplygin there is also a coordinate system tp, y where 
energy has the following form 

H = i(y) {p% + p 2 y ) + 0{y) cos tp 

for some functions 7 and /3, see 0. Since H = j{y){p 2 + p 2 ) is the Hamiltonian of 
the geodesic flow of a smooth metric on S 2 , this coordinate system is unique up to 
linear combinations. 

On the other hand there is an integral F = p^+np^H+G^p^^Py, <p, y) where constant 
k^O and G(p v ,p y , ip, y) is a polynomial linear in momenta (A.V. Bolsinov, private 
communication), see also 0. 

If our series includes the case of Goryachev-Chaplygin then it follows that in the 
case of Gorychev-Chaplygin there are two independent first integrals. So, our series 
of integrable conservative systems on S 2 does not include the case of Goryachev- 
Chaplygin. □ 

Commentar. Numerically we got that the value of T is equal to ~0, 57735. 
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